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We make use of a semiclassical method for calculating the suppression exponent for topology 
changing transitions in high-energy electroweak collisions. In the Standard Model these processes 
are accompanied by violation of baryon and lepton number. By using a suitable computational 
technique we obtain results for s-wave scattering in a large region of initial data. Our results show 
that baryon and lepton number violation remains exponentially suppressed up to very high energies 
of at least 30 sphaleron masses (250 TeV). We also conclude that the known analytic approaches 
inferred from low energy expansion provide reasonably good approximations up to the sphaleron 
energy (8 TeV) only. 
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I. INTRODUCTION 

Non-perturbative phenomena related to tunneling are 
often encountered in quantum field theory. Well known 
examples are false vacuum decay and instanton-like tran- 
sitions (the latter are accompanied by non-conservation 
of fermion quantum numbers). When these phenomena 
are governed by a small coupling constant they can gen- 
erally be studied by semiclassical methods. This is cer- 
tainly the case at low energy or in situations which in- 
volve large number of quanta in the initial state. At 
low energy, collision processes can be well described by 
a semiclassical approximation relying on the existence of 
classical Euclidean time solutions to the equations of mo- 
tion interpolating between initial and final states. In the 
examples mentioned above these are bounce and in- 
stanton Q , respectively. The probability of the process 
is then proportional to the exponent of the Euclidean 
action of these solutions. As the action is inversely pro- 
portional to the (small) coupling constant, the processes 
are highly suppressed. The effect of low energy excita- 
tions in the initial state (colliding particles) gives only a 
pre-exponential factor and is inessential. 

The situation changes at high energy, namely, at en- 
ergy of the order of the tunneling barrier height which 
separates initial and final states. In general, there ex- 
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ists a static unstable solution to the equations of mo- 
tion that lies on top of the potential barrier Q (properly 
speaking, at a saddle point of the potential) . In field the- 
ory this solution is often referred to as the "sphaleron" , 
a name which we will use throughout this paper. The 
minimum height of the barrier is precisely the sphaleron 
energy Esph- Naively, from the analogy with quantum 
mechanics of one degree of freedom, one would expect 
that at energy higher than the sphaleron energy the ex- 
ponential suppression disappears. This is indeed what 
happens at finite temperature [l, |p , | p. B, H, M, ITol Ull IT^ ■ 
finite fermion density 0, IT^ llsL Il6l , or in the pres- 
ence of heavy fermions in the initial state [l^ . IT^ l20l | . 
But in high energy particle collisions this is not neces- 
sarily the case, due to the fact that the characteristic 
size of the sphaleron configuration is much larger than 
the wavelength of the incoming particles. At the same 
time the application of a semiclassical technique becomes 
problematic because the initial state no longer involves a 
large number of quanta. 

As was first noted in Refs. [2ll [2^. at relatively low 
energy the corrections to the collision-induced tunnel- 
ing rate can be calculated by perturbative expansion in 
the background of the instanton. Further studies showed 
that the actual expansion parameter in most models, in- 
cluding electroweak theory, is E/Esph pL I2I I2ll2^ 
and the total cross section of the induced tunneling has 
an exponential form 



atot{E) ~ exp —FHGiE/E,ph) 
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where g is the small coupling constant and the function-"^ 
Fhg{E / EspY^) is a series in fractional powers of E/Esph 
(for a review see [IS S 113). 

While the perturbation theory in E/E^ph is limited to 
small E, the general exponential form of the total cross 
section implies that there might exist a semiclassical-typc 
procedure which would allow, at least in principle, to cal- 
culate Fhg{E / Esph) a.t E > Egph- However, since the 
initial state of two highly energetic particles is not semi- 
classical, the standard semiclassical procedure does not 
apply and a suitab le g eneralization is needed, which was 
proposed in Refs. [soi |^ IS3 ^^^"^ further developed in 
Refs. [s^jHI- The corresponding formalism reduces the 
calculation of the exponential suppression factor to a cer- 
tain classical boundary value problem, whose analytical 
solution is not usually possible. 

This approach is based on the conjecture that, with 
exponential accuracy, the two-particle initial state can 
be substituted by a multiparticlc one provided that the 
number of particles is not parametrically large (although 
not proven rigorously, this conjecture was checked in sev- 
eral orders of perturbation theory in E/Egph in gauge 
theory [sil Issj and explicitly in quantum mechanics with 
two degrees of freedom [Sa, |33)- The few-particle ini- 
tial state, in turn, can be considered as a limiting case 
of a truly multiparticle one with the number of particles 
N = N ■ {Air/g'^), when the parameter N is sent to zero. 
For the multiparticle initial state the transition rate is 
explicitly semiclassical and has the form 

aiE,N)^expi^~^F{E/E,pi,,N)Y (1) 

According to the above conjecture, the function 
Fhg{E I Esph), corresponding to the two-particle incom- 
ing state, is reproduced in the limit TV ^ 0, 

lim F{E/Esph,N) = FHoiE/Esph). 

Therefore, albeit indirectly, the function Fhg{E / Esph) 
is also calculable semiclassically. 

Within the semiclassical framework, the function 
F{E/Esph,N) is determined by the action evaluated at 
a particular solution to the classical field equations [s^ 
on a certain contour in complex time. In this formula- 
tion, the problem, at least in principle, is amenable to a 
computational solution. Namely, one has to solve the cor- 
responding classical boundary value problem numerically 
and calculate the function F{E / Esph, N), which then can 
be used to extract information about Fhg{E/ Esph)- 

The implementation of this technique is nevertheless 
highly non-trivial. The differential equations one en- 
counters are partially of the hyperbolic type (along the 
Minkowskian parts of the time contour) and partially 
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of elliptic type (along the Euclidean part), which makes 
their numerical solution particularly challenging. In the 
electroweak theory, additional difficulties arise from the 
need to deal with the large number of internal degrees of 
freedom, unphysical modes due to gauge invariance, and 
time translational symmetry which cause an unwelcome 
degeneracy in the numerical procedure used to find the 
semiclassical solutions. 

Moreover, at high energy (roughly, at energy higher 
than the sphaleron energy) tunneling solutions, interpo- 
lating between vicinities of different vacua in finite time, 
cease to exist. This subtlety turned out to be a general 
problem in the description of tunneling in systems with 
many degrees of freedom, and it has to do with the non- 
trivial way tunneling occurs at high energy — the system 
prefers to create a state close to the sphaleron, which 
then decays into the correct vacuum. To find the cor- 
responding suppression exponent numerically one has to 
use a properly regularized version of the boundary value 
problem, developed in Ref. [38l |. 

In a long program of investigations we have been able 
to gradually overcome all of these hurdles. Preliminary 
results for energies below the sphaleron energy were re- 
ported in Ref. [33j . Now we are in a position of presenting 
what we are confident is the full solution to the numerical 
problem for a wide range of energy and incoming parti- 
cle number, including energies above the sphaleron. The 
field configurations we analyse in this paper are restricted 
to be spherically symmetric in space. Hence our results 
apply, strictly speaking, to s-wave scattering only. 

In this paper we will concentrate on obtaining the 
suppression exponent for collision-induced tunneling in 
SU{2) gauge model with the Higgs mechanism, corre- 
sponding to the electroweak sector of the Standard Model 
at zero 9w- The problem is particularly interesting be- 
cause of the baryon number violation which accompanies 
such processes |39| and the relatively low sphaleron en- 
ergy Esph — 8TeV. Though computational limitations 
do not allow to reach literally zero value of the rescaled 
number of particles N, corresponding to particle colli- 
sions, we were able to extrapolate the results to zero N 
and get a bound on the suppression exponent (strictly 
speaking, for s-wave scattering) and also provide an esti- 
mate for this exponent. 

In sectioninwe present the detailed formulation of the 
problem, outline the method and present the main phys- 
ical results. In section iHll we give the derivation of the 
semiclassical method for the gauge model. The lattice 
formulation of the equations and subtleties appearing in 
the discretized version are given in section IIVI Applica- 
tion of the regularization method of Ref. [Sg is described 
in section Detailed numerical results are presented in 
section IVTl Our conclusions arc in section IVIII 
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II. FORMULATION OF THE PROBLEM AND 
MAIN RESULTS 

Non-Abelian gauge models have an infinite number of 
topologically distinct vacua, labeled by an integer topo- 
logical number. Processes changing the topological num- 
ber arc accompanied by violation of fcrmion (baryon and 
lepton) numbers |39| , a phenomenon of great interest for 
cosmology and particle physics. The topologically dis- 
tinct vacua are separated by a potential barrier, whose 
height, in models with the Higgs mechanism, is given by 
the sphaleron energy. Topology changing transition may 
occur via tunneling at low energies or, at sufficiently high 
energy and suitable initial state, via classical evolution 
over the sphaleron. 

In this paper we study a four-dimensional model which 
captures all the important features of the Standard 
Model — an SU{2) gauge theory with the Higgs doublet. 
This model corresponds to the bosonic sector of the Stan- 
dard Model with 9w ~ 0. To the leading order in the 
coupling constant, the effect of fcrmions on the gauge and 
Higgs fields dynamics can be ignored 0| . The action of 
the model is 



S 



4™./ 1 2 



(2) 



A(*1'$ - l)^! , 



where 



dp. An dyAp 

{dp - lAp)^ 



i\Ap-, Ajy] 



(3) 
(4) 



p — A^a°-/2 and aw = ff^/47r. Here we have 



with A 

eliminated inessential constants by an appropriate choice 
of units. The dimensional parameters can be restored 
noting that in the normalization the gauge boson 
mass is 



w 



1 

V2 



(5) 



and the Higgs boson mass is 



M 



H 



8AMw 



In most of our calculations the Higgs self-coupling A was 
set equal to A = 0.125, which corresponds to Mh = Mw- 
The dependence on the Higgs boson mass is very weak, 
so this is a reasonable approximation. Also we often omit 
the omnipresent overall factor 1/aw- 

Any vacuum configuration in this model can be ob- 
tained from the trivial vacuum A^ = 0, $ = <f>t, = ^ 

by a certain gauge transformation U{x). We will be us- 
ing the temporal gauge = 0, where the vacuum con- 
figurations are described by time- independent U{x), cor- 
responding to residual gauge invariancc. In this gauge 



field values at spatial infinity cannot change during the 
evolution (otherwise the kinetic term becomes infinite) 
and thus one considers only those U {x) which have some 
fixed asymptotics at spatial infinity. Often the asymp- 
totic U{x — > oo) — > 1 is used, so any vacuum configura- 
tion corresponds to a mapping from space with iden- 
tified infinity, which is homotopically equivalent to S'^ , to 
the gauge group SU{2) ^ . The degree of this mapping 
is precisely the topological number of the corresponding 
vacuum. A gauge choice of this form is convenient for 
analysis of the excitations about the trivial vacuum. For 
other purposes it may however be useful to choose an 
alternative behavior of the gauge function at spatial in- 
finity, like U{x) — > exp{jcra;/|a;|}, which maps the S"^ of 
spatial infinity to the equatorial of the SU{2). The 
two neighboring vacua then map the space either to 
north or south hemisphere of the SU{2). In this gauge, 
the sphaleron configuration takes the simplest form, and 
we will use this gauge everywhere in this paper, except 
for the analysis of the mode expansion in the initial state. 

Numerous perturbative attempts were made to find 
the probability of the collision-induced topo logy chang- 
ing transitions in this model (see Refs. [2^ I2M l29l | for 
reviews), giving reliable results only for relatively low 
energies. A non-perturbative study of classically allowed 
over-barrier transitions was presented in Ref. |4ll |. All 
solutions found in Ref. |4lj | are configurations with large 
numbers of particles in the initial state and thus they do 
not correspond to realistic collisions. Another approach, 
pursued in this pap er, is to use the semiclassical method 
of Refs. jsO, I31I I32L Isl l adapted for theories with gauge 
degrees of freedom. This method was implemented in 
Ref. [33, where the results were obtained for energies 
below i?sph, what suggested that at the sphaleron energy 
the suppression is still strong. However, a straightfor- 
ward application of the technique of Refs. [30II311 . 13^ Is^ 
fails for energy above the sphaleron due to the problems 
one encounters as the energy approaches the height of the 
barrier in systems with many degrees of freedom. These 
problems were studied in detail in the context of a quan- 
tum mechanical model in Ref. ;38<] , where a regularization 
technique was suggested to overcome them. 

The basic idea in the proposal of Refs. [30ll3lll3^ l3^ 
is that, instead of a process with exclusive, two-particle 
initial state, one considers a topology changing process 
with inclusive initial state characterized by definite en- 
ergy E and incoming particle number N. The transi- 
tion probability (t{E, N) can then be used to provide a 
bound on the exclusive two-particle cross-section, while 
the two-particle transition exponent is obtained in the 
limit awN — > 0. 

The inclusive probability of tunneling from a state with 
fixed energy and number of particles is 



a{E,N)=J2\(f\^PEpN\^) 

ij 



(6) 



where S is the 5'-matrix, Pe^n are projectors onto sub- 
spaces of fixed energy E and fixed number of particles 
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restoring it only in the final results. The boundary term 
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FIG. 1: The contour in complex time plane used in the for- 
mulation of the boundary value problem @. Crossed circles 
represent singularities of the field. If the field is spherically 
symmetric in space, the singularities closest to imaginary axis 
occur at r = 0, for other r the singularities generally move to 
larger | Ref|. 



N, and the states \i) and |/) are perturbative excitations 
about topologically distinct vacua. This matrix element 
can be written in double path integral representation. 
For large N = N /aw and E = E/aw the path integral 
can be calculated in the semiclassical approximation, and 
this leads to the problem of solving the equations of mo- 
tion of the system on a special contour in complex time 
plane, which detours around singularities, as shown in 
Fig. The presence of branch cut singularities can be 
inferred from the following argument. One notices that, 
for energy below the sphaleron energy, if one continues 
the solution along a line parallel to the real axis, be this 
via a forward integration of the equation of motion from 
the AB part of the contour or a backward integration 
from the CD part of the contour, the field must fall back 
to the original topological sector. On the other hand, by 
construction, on the AB and CD parts of the contour the 
solution must be in different topological sectors. Thus 
the solution must also be in different topological sectors 
on the AB part of the contour and on the negative real 
axis and, likewise, on the positive real axis and the con- 
tinuation of the AB segment to positive time. This may 
happen only if two branch cut singularities exist on the 
two sides of the BC part of the contour (see Fig.P). 

Eventually the semiclassical approximation produces 
the following result ((^ here stands for all physical fields 
in the model). 



Ait 

(j{E, N) - exp <j F{E, N) 

aw 



(7) 



47r ~ ~ 

F{E, N) = 2 Im Sabcd{v) -N9-ET-Re B, 

aw 



Here Sabcd{'^) is the action along the time contour, the 
parameters T and 9 are Legendre conjugate to E and 
N; the parameter T is the same as in Fig. ^ we will 
have to say more about 9 later on. In what follows we 
will usually drop the tilde over the rescaled energy and 
incoming particle number, and the overall 1/aw factor. 



6,; = ; /dfe(.4/_fce-2--(^'-^^/2)-5*g::^e2-'"(^'-'^/2)) 



is written using frequency components fj^ and g^. of the 
field on the part A of the contour: 



dk 



(8) 



-iu)k{t-iT /2)+ikx 



The field <p satisfies the field equation 

— -0 

5if 



(9a) 



At initial time the frequency components of the solution 
should satisfy the following equation ("0 boundary con- 
dition" ) 



fk = e 5fc 



(9b) 



For 9 different from zero this equation implies that the 
field must be continued to complex values. For a complex 
field, like <I> in its real and imaginary parts must be 
continued to complex values separately. 

On the final part of the contour (CD), the field must 
satisfy the reality condition 

lm(p{x,Tf ^ oo) , lm(p{x, Tf oo) ^ (9c) 

(for complex fields, such as <t> in (j2J), this means that both 
($ + $t)/2 and ($ - $^)/2i must be real). 

Equations (|9a|l - l|^ specify the boundary value prob- 
lem corresponding to the induced topological transition. 

The equations obtained by variation over the auxiliary 
parameters T and 9 are 



E^ / dkujf^f^gl 



N 



(10) 

(11) 



These equations indirectly fix values of T and 9 for given 
energy and number of particles. Alternatively, one can fix 
T and 0, solve the boundary value problem ^ and obtain 
the corresponding values of E and N using and 
This is especially convenient in numerical calculations. 

The interpretation of the solutions to the boundary 
value problem Q is as follows. On the part CD of the 
contour, the saddle-point field is real asymptotically; it 
describes the evolution of the system after tunneling. On 
the contrary, it follows from boundary conditions (|9bl) 
that in the initial asymptotic region the saddle-point field 
is complex whenever 9^0. Thus, the initial state which 
maximizes the probability Q is not described by a real 
classical field, i.e. this stage of the evolution is essentially 
quantum even at TV ~ 1/aw- 
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FIG. 2: Regions in the E-N plane. 



There is a subtle point concerning the boundary condi- 
tion It can be satisfied in two different ways. Either 
the solution is exactly real on the whole CD part of the 
contour and is close to vacuum, or it has an exponentially 
decaying imaginary part and approaches the sphaleron 
along the complexified unstable direction. This subtlety 
is important for the analysis at high energies fsection lV)! . 

The solutions to the boundary value problem can be 
found numerically for different values of E and N . In this 
paper we study solutions that have spherical symmetry 
in space. One expects that these are most important 
for large enough N\ perturbative calculations about the 
instanton suggest that spatial spherical symmetry is rele- 
vant at relatively low energies and all N . We do not have 
a convincing argument in favor of spherical symmetry for 
few particle collisions at very high energies; in any case, 
our results as they stand, are valid for s-wave scattering. 

Our numerical analysis shows that the E — N plane is 
divided into several different regions (see Fig.[2Jl. Values 
of i? < • min(Mvi/, Mh) are trivially excluded by kine- 
matics. For relatively low energies (region A) the transi- 
tions between the topologically distinct vacua can occur 
only via tunneling. At the sphaleron energy E'sph the 
situation changes. A slight excitation of the sphaleron 
along the unstable direction gives origin to a solution of 
the classical equations of motion which evolves towards 
different topological sectors at large negative and posi- 
tive times. Since the sphaleron has exactly one negative 
mode, there is only one infinitesimal deformation of this 
type, and thus the corresponding solution has definite 
number of particles A'sph in the initial state. At higher 
energies one may add excitations of the positive modes 
above the sphaleron to obtain over-barrier solutions with 
different, and, in particular, smaller initial particle num- 
ber. These solutions belong to the domain of classically 
allowed transitions (region B in Fig. ^ , where the topol- 
ogy changing processes are unsuppressed. The boundary 
between region A and region B corresponds to configu- 



rations staying for an infinite time close to the sphaleron 
(since there are no bound states in the sphaleron back- 
ground 1 £^11 excitations about the sphaleron fly away 
at finite time, so the field relaxes to the sphaleron solu- 
tion). 

In the classically forbidden region A there is a spe- 
cial family of solutions, corresponding to = in the 
boundary value problem. These are represented by the 
line Epi{N) in Fig. |2 In this case, the boundary con- 
dition IjQbp reduces to the reality condition imposed at 
Imt — T/2. The solution to the resulting boundary value 
problem is the periodic instanton of Ref . ^3 ■ The peri- 
odic instanton is a real periodic solution to the Euclidean 
field equations with period T and turning points at t = 
and t = iT/2 (modT). When analytically continued in 
the Minkowskian direction through the turning points, 
the periodic instanton is real on the lines Imt = and 
Imt = T/2 and therefore satisfies the boundary value 
problem Q with 6 — 0. Like any other solution lineariz- 
ing at large negative times at part A of the contour of 
Fig. ^ the periodic instanton has a certain number of 
incoming particles, Eq. For given energy E below 

the sphaleron, this number is such that the suppression 
exponent F{E,N) has a minimum, i.e., the transition 
occurs at maximum rate. 

The classically forbidden region A is further subdivided 
into two regions. For low energies (region A.I) the system 
is close to the vacuum on the final part of the evolution, 
so the boundary condition leads to the exact reality 
of the fields on the part CD of the time contour. At en- 
ergies higher than the sphaleron energy (precisely, on the 
right of the line Ei{N)) the system ends up close to the 
sphaleron (with extra outgoing waves in the sphaleron 
background). In this case Eq. is truly asymptotic. 
So, the system tunnels "on top" of the barrier, creating 
an unstable sphaleron configuration, which then decays 
with probability of order 1 to any of the two neighbor- 
ing vacua. This situation is realized in the region A.II. 
This new qualitative feature of the tunneling at high en- 
ergies emerges from the existence of the bifurcation of 
the solutions and is not seen in any order of perturbative 
expansion around the instanton. Non-perturbative ap- 
proaches, however, capture this feature (see also Ref. |4J| 
for similar results in the context of false vacuum decay). 

Our numerical results for the suppression exponent 
in the whole classically forbidden region are presented 
in Figs. 13 0] The almost vertical line in Fig. O sepa- 
rates the two regions (denoted by A.I and A.II in our 
earlier discussion) where the tunneling process assumes 
characteristically distinct features. It also represents 
the frontier beyond which numerical calculations based 
on a straightforward implementation of the method of 
Rcfs. j30j 31, 32, 34] appear to fail. It is clear from Fig.|2| 
that our improved numerical technique can go well be- 
yond that frontier. Reference [45| presents a comparison 
between our results and the analytic predictions for the 
suppression exponent F{E^ N) in the limit of small en- 
ergy. The two are in remarkable agreement which pro- 
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FIG. 3: Lines of F(E, N) — const. Lines are labeled by 
the values of the suppression exponent —aw log a = AttF. 
Diagonal line directed from the sphaleron towards the origin 
is the line of periodic instantons. Energy E is in units of 
Mw/cew, number of particles A'' is in units of 1/aw- The 
line labeled hy {F = 0) is the boundary of the classically 
allowed region E = Eq{N). The "fuzzy" line represents the 
approximate boundary of the classically allowed region found 
in over-barrier calculations of Ref. |4l|. 
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FIG. 5: Lower bound on the suppression exponent for two- 
particle collisions, dashed and dashed-dotted lines. Dotted 
line is the estimate of Refs. Ea, ES. 
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FIG. 6: Estimate of the suppression exponent for two-particle 
collisions Fhg{E) (solid line), lower bound on Fhg{E) 
(dashed line), low energy analytic prediction (1121 (rare dotted 
line) and analytic estimate of Refs. |4^. li^ (dotted line). 



FIG. 4: Dependence of the suppression exponent on the num- 
ber of particles A'' for different energies. Numbers near the 
curves are energies in units of Mw/ctw- 



vides a gratifying check of the numerical calculations. 

Another interesting comparison can be made with the 
results of Ref. where the real-time overbarrier so- 
lutions close to the boundary of the classically allowed 
region were searched via Monte-Carlo techniques^. In 



this way, an approximation (which, at the same time, is 
an upper bound) for the boundary of the classically al- 
lowed region was obtained. It is seen that the results of 
Ref. ^3 are reasonably close to the boundary Eq{N), 
found in our calculations. 

Our results by themselves do not reach the physically 
interesting N = limit (corresponding to particle colli- 
sions) . Studying lower N in numerical calculations would 
need lattices with larger number of lattice points and 
would require quite substantial amounts of time even on 



^ In Ref. |4ll the coupling constant A was chosen 0.1, while we 
use A = 0.125. We performed a set of calculations for A = 
0.1. The dependence on A is so weak, that the difference for the 
results would be invisible in the graph. Much larger discrepancies 



appear because of the diff erent lattice parameters used in the two 
calculations. (In Ref. |41| . having only to solve for the real time 
evolution of the fields, it was possible to use a larger lattice and 
a finer lattice spacing than in the present calculations.) 
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powerful present day supercomputers. Therefore, some 
extrapolation must still be used to get insight on the 
suppression factor for actual particle collisions. As we 
seek such extrapolations, we notice first that it is quite 
straightforward to obtain a lower bound on the suppres- 
sion exponent F. Insofar as 9 increases as iV ^ 0, and 
{ATi)dF / dN = —9, by simply continuing F with a linear 
function of A'^ for each energy one obtains a lower bound 
on F. This bound is shown in Figs. HJ dashed line. It 
indicates that up to the energy SMiy/aw — 20TeV the 
suppression is still high: the suppression factor is smaller 
than e-6° - lO'^^ for aw ~ 1/30. 

For very high energies a bound may be constructed 
by exploiting the observation that the lines of constant 
F in E — N plane have positive curvature (see Fig. ISJ. 
So, by extrapolating these lines linearly to TV = one 
obtains another lower bound on the suppression exponent 
F{E,N = 0). This bound is displayed in Fig.|31 dashed- 
dotted line. One can see that exponential suppression 
continues up to an energy of at least 250TeV. 

One may also attempt to estimate the function F{E) 
itself. As we discuss in Section IVII a good estimate is 
obtained by extrapolating, instead of F{E,N), the func- 
tion T{N) at fixed energy, as T{N) is approximately 
linear in A^. Up to the sphaleron energy, the estimate 
obtained in this way is close to the one loop analytic 
result ji^ \^ Is^, 1^, which gives three terms in the 
low-energy expansion. 



47r , , 47r 

F(E) = — 

aw ctw 



9fE_ 

8 \Eo 



4/3 



9 / E 



16 \En 



(12) 



where Eq = ^/6nMw /ctw ■ Below the sphaleron, our 
estimate is also consistent with the analytic estimate 
of Refs. |46l l47j . On the other hand, the behavior of 
Fhg{E) changes dramatically a.t E > E^ph- We at- 
tribute this to the change in the tunneling behavior — 
at E > Esph the system tunnels "on top of the bar- 
rier" . Our numerical data show that the suppression ex- 
ponent Fhg{E) flattens out, and topology changing pro- 
cesses are in fact much heavier suppressed at > -Bsph 
as compared to the estimate ()12|l and the estimate of 
Rcfs. We show our estimate, together with ana- 

lytical estimates and our lower bound, in Fig. |S1 

It is worth noting that similar effects of dramatic 
change of the behavior of the system at high energies 
were observed in lattice calculations of instanton distri- 



bution in QCD in |5al53|. 

Thus, our numerical results, albeit covering a limited 
range of energies and initial particle numbers, enable us 
to obtain both lower bound for and actual estimate of 
the suppression exponent for the topology changing two- 
particle cross-section in the electroweak theory well above 
the sphaleron energy. This cross section remains expo- 
nentially suppressed up to very high energies of at least 
250TeV. In fact, the energy, if any, at which the exponen- 
tial suppression disappears, is most likely much higher, 
as suggested by comparison of our lower bound and ac- 
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FIG. 7: The contour used to derive the boundary value prob- 
lem. 



tual estimate at energies exceeding significantly -Egph, see 
Fig.El 



III. THE METHOD OF RST 

A. General formulation 

The quantity we wish to calculate is (t{E, N), the prob- 
ability of transition from a state with fixed energy E and 
number of particles N about one vacuum to any state 
about another vacuum. The method of semiclassical cal- 
culation of this inclusive multiparticle probability was 
formulated in Refs. (aSEHHlHl. We call it the RST 
method for brevity, and here we review this prescription 
in brief. 

The inclusive multiparticle probability ((BJ can be writ- 
ten in functional integral form, where the semiclassi- 
cal approximation is equivalent to the saddle-point in- 
tegration. The double path integral representation for 
a{E, N) reads [13 



a{E,N) = J dMTda^daldbkdbld(p{x)d(p' (x) exp^^-iNd - iET - J dka^c 



le-'^-'^"'^ - I dkbkbl 



+ B,{ak.^,) + Bfibl, Vf) + B*{a*_k, + i?;(6_fc, ^'^) + iS{^) - iS{^')] . (13) 
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Here Lp stands for all physical fields of the theory. The boundary terms Bi and Bf are 

1 



dk 



1 



(14) 



r 



where ipi jik) are the spatial Fourier transforms of the 
field at initial and final times Ti and Tj, respectively. 
The limit Tij — > =poo is assumed at the end of the cal- 
culation. The complex integration variables and 
come from the coherent state representation of initial and 
final states; they are the classical counterparts of annihi- 
lation and creation operators. The integration over these 
variables implements the summation over initial and final 
states in (0. The functional integrals over 1^9(2;) and ^p'ix) 
come from the amplitude and complex conjugate ampli- 
tude, respectively. The integrations include the bound- 
ary values Lpi f and Lp'^ ^. Integration over T and 1? serve 

I 



to project onto the subspaces of fixed E and N, respec- 
tively. 

The integral (U^ can be evaluated in the saddle point 
approximation, as long as the exponent is proportional to 
l/aw, implicitly present in the expression, and N,E 
1/aw- 

Let us now discuss the saddle-point equations for the 
integral We will see that these equations reduce to 
a certain boundary value problem for the fields ip and 
ip' . The variables a^, a^, 6^ and 6^ enter the exponent 
quadratically and can be integrated out, yielding 



a{E,N) = J dMTdLp{x)dLp'{x)'[lS{ipf{k) - ip'f{k)) x expl^~iNd ~ iET + iS{Lp) - iS{ip') 

2-((l+7D[¥'«(fc)^.(-fe) + ^:(fc)^U-fc)]-47fc^.(fc)^-(-fc))| , (15) 



2 ■ ''"-^ 



r 



where 



Ik 



An important feature of the representation H15(l is that 
the exponent in the r.h.s. contains only the action and 
the boundary terms. Thus, the discretization of this ex- 
ponent is relatively straightforward. 

Let us turn to the saddle point equations. Varying the 
exponent with respect to the fields ip{x) and ^'[x) we 
find 



Sip 5ip' ' 



(16) 



i.e. the usual field equations. The boundary conditions 
for these equations come from the variation with respect 
to the boundary values of the fields. At f = Tj, be- 
cause of the (5-function, the variations are subject to 
the constraint 5iff{x) = 5ip'j{x) (at T/ 00). Since 
SS/6(p{Tf,x) = Lp{Tf,x) we obtain 



^{Tf,x) - ip'{Tf,x) 



(17) 



Thus, in the final asymptotic region the saddle-point 
fields ip and </?' coincide. 



The variation with respect to ipi and ip[ leads to two 
equations which can be written in the following form, 

i(pi{k) + uJk^Piik) = (i^iik) + uJk^P^ik)) , 

-iipi{k) + uJkipi{k) = — (-iip'^{k) + LUkiPiik)) . 
Ik 

These initial boundary conditions simplify when written 
in terms of frequency components. In the initial asymp- 
totic region (t — *■ —00), where ip and ip' are free fields, we 
can write 



ip{x) = 



dk 

-T- 



■.{fk 



-iLOtst'\-ikx 



* iuji^t—ikx 



^p'[x) = / ^ -^u>^t+^kx 

Then the initial boundary conditions become 



}, (18) 



iuji^t — ikx 1 



fk 
9k 



Ikf'k 



Ik 



I* 
-9 k 



(19) 



Finally, there are two saddle-point equations which come 
from the variation of the exponent in (|15|l with respect 
to § and T. These equations determine the saddle-point 
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values of and T as functions of E and N. In terms of 
frequency components fk and gp, they read (after using 
boundary conditions ) 

E ^ J dkuj^gl, (20) 

N = j dkf^gl (21) 

One may recognize the usual expressions for the energy 
and the number of particles contained in the free classical 
field, Hf, = ff,gl. being the occupation number in the 
mode with spatial momentum k. 

The field {x) originates from the complex conjugate 
amplitude. This suggests that its saddle point value is 
complex conjugate to that of ^{x). Indeed, the Ansatz 

is compatible with the boundary value problem H16|l - 
(|21|l . Then the saddle point values of T and d arc pure 
imaginary 

T = iT , i9^ie , 

provided the initial energy H20|l and particle number H21|l 
are real. The boundary conditions H17|l imply then that 
the field ip is real asymptotically at final time 

lmip{Tf,x) ^ , lmLp{Tf,x) ^ for T/ ^ +oo , 

while Eq. H19() relates the positive and negative frequency 
components of the field ip in the initial asymptotic region 

fk = 7fc5fc > 

where 

7fc - 0-"-^-^ , (22) 

Until now, the initial time Ti was real. However, it is 
convenient to reformulate the boundary value problem 
directly in terms of the fields on the contour ABCD, at 
which the initial time has imaginary part ImT"; = T'/2 
(see Fig. [7||. The analytical continuation in the initial 
asymptotic region can be done explicitly by means of 
Eqs. ((TH|) . In Eqs. lfT^ - l(^ this continuation results in 
the substitution of 7^ by 

7fc = e '"^ '■ 

The simplest boundary conditions are obtained in the 
case when the contour height in imaginary time T' is 
equal to the parameter T, leading to fc-indepcndent 7 



boundary value problem (|16() " (|21|l is equivalent'^ to ©- 
(|ll|l . This is the boundary value problem we solve nu- 
merically in the present paper. 

Let us discuss some subtle points of this boundary 
value problem. First, one notices that the condition of 
asymptotic reality does not always coincide with the 
condition of reality at finite time. Of course, if the solu- 
tion approaches the vacuum on the part CD of the con- 
tour, the asymptotic reality condition (|9c|l implies that 
the solution is real at any finite positive t. Indeed, at 
large enough time the system evolves in that case in the 
linear regime, so the condition H9c|) means that all phys- 
ical modes should be real. Due to the equations of mo- 
tion the fields are then real on the entire CD-part of 
the contour. This situation corresponds to the transition 
directly to the neighboring vacuum. However, the situ- 
ation can be drastically different if the solution on the 
final part of the time contour remains in the interaction 
region, i.e. close to the sphaleron. Since one of the ex- 
citations about the sphaleron is unstable, there may ex- 
ist solutions which approach the sphaleron exponentially 
along the complexified unstable direction. In that case 
the solution may be complex at any finite time, and be- 
come real only asymptotically, as t — > -|-oo. Such solution 
corresponds to tunneling to the sphaleron; afterwards the 
system rolls down classically to the correct vacuum (with 
probability of order 1, inessential for the tunneling expo- 
nent F). We will see in section |V] that the situation of 
this sort indeed takes place at high energies E > £^sph- 

Second, the initial boundary conditions H19|l (imposed 
on the real time axis) mean, that (p and </?' ~ (p* are 
different at large negative time, while at large positive 
time they coincide because of the condition (|17|l . For 
solutions ending in the vacuum at positive time (so that 
the fields are exactly real at finite t > 0), this means that 
there should exist a branch point in the complex time 
plane: the contour in Fig. ^ winds around this point and 
cannot be deformed to the real time axis. This argument 
does not work for solutions ending on the sphaleron at 
t — > -l-oo, so branch points between the AB-part of the 
contour and the real time axis may be absent. We have 
found that this is indeed the case at high energies (cf. 
Ref. [H). 



B. Reduction to spherically symmetric 
configurations 

Here we consider spherically symmetric configura- 
tions [S^l of the SU (2)-Higgs theory. The reason is that 
one can entertain the expectation that the most impor- 
tant tunneling configurations possess maximum spatial 
symmetry. On the other hand, without the simplifica- 



In this case one arrives at the boundary condition (|9b|) 
and the contour ABCD with height T/2 shown in Fig.^ 
This formulation will be used in most cases. Then the 



^ The boundary term in is obtained from the boundary terms 
in I15i by making use of the 6 boundary conditions 19bl . 
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tion provided by spherical symmetry the computational 
cost of the numerical analysis would be prohibitive. 

In the spherically symmetric Ansatz the original fields 
are expressed in terms of six real two-dimensional fields 
ao, ai, a, /3, fj, and v as follows 



Aoix,t) 



1 . N 

-ao(r, t)(T • n 



ai(r, t)cr • nrii H ^ — (ai — cr • rtrii) 



l + /3(r,t) 



(23) 



$(a;, t) = [/i(r, i) + iiy{r, t)a- ■ n]^ , 



where n is the unit three-vector in the radial direction 
and ^ is an arbitrary constant two-component complex 
unit column. This Ansatz is symmetric under spatial 
rotations complemented by appropriate rotations in the 
gauge group and custodial global symmetry transforma- 
tions. The action |5J) expressed in terms of the new fields 
becomes 



S 



dt I dr 



I? 



(24) 



r 



where the indices /i, v run from to 1 and 



The equations of motion obtained from H24|l are 



X = a + if3 x = a-«/3 

(f> ~ II + (j) ~ II — iv 

D^X= {dt,-iaf,)x D^x = [dt, + ia^)x 



(25a) 
(25b) 
(25c) 
(25d) 

(25e) 



Note that the overbar on (j>, x and -D^ denotes changing 
I — > — i in the definitions (|25|l above, which is the same 
as complex conjugation only if the six fields a^, a, /?, /j, 
and V are real. In the boundary value problem Q these 
fields become complex and overbar no longer corresponds 
to normal complex conjugation. 



(26a) 



9o(rVoi) - *[X^1X - X^iX] + ^r^[0i^ - 4>Di 



D^r'D^ + -ixx+l) + ^^r 
D^r'D^ + l{xx + l) + '2Xr' 









(26b) 


l) + i<^0 


X = - 


2^ 


(26c) 


1) + i# 


X = - 


2^ 


(26d) 


+ 2Xr^{4><, 


6-1) 


(j) = ix(f> 


(26e) 



1) 



(f) = iX(t> 



(26f) 



Equation l|26a(l is of the first order in time — it is Gauss' 
law. 

The spherical Ansatz (I23II has a residual U{\) gauge 
invariance 



X ^ c^^'x 



(27a) 
(27b) 
(27c) 



with gauge function VL{r,t). The complex "scalar" fields 
X and (j) have U{1) charges 1 and 1/2 respectively, 
is the C/(l) gauge field, /^^ is the field strength tensor, 
and in (|25|l is the covariant derivative. The residual 
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(7(1) gauge invariance must be fixed when solving the 
equations numerically. We choose the temporal gauge 
ao = 0. In this gauge, if Gauss' law is obeyed at some 
moment of time, the other five equations guarantee that 
it is obeyed at any time. This means, in fact, that one 
of the equations is redundant, and one of the fields is 
not physical — it can be expressed in terms of the other 
four fields and their derivatives using Gauss' law. How- 
ever, numerically it is easier to solve five second order 
equations of motion imposing Gauss' law as one of the 
boundary conditions. Also, in the oo = gauge, there 
remains a gauge freedom with time independent gauge 
function, and this should also be fixed by boundary con- 
ditions. 

The trivial space-independent vacuum of the model is 



r = oo. 



±1 







(28) 



Other vacua are obtained from the trivial one by the 
gauge transformations 



a 



fi vac 9^ 



in 



Xvac 

/ J in/ 2 

0va,r = ±e ' 



(29a) 
(29b) 
(29c) 



By regularity, should be zero at the origin. Vacua 
with different winding numbers correspond to — > 2mr 
as r ^ oo. For such values of fi, the fields of the original 
four-dimensional model are constant at spatial infinity, 
and this is the standard choice. It allows for a simple 
description of the topological properties of vacua: since 
the sphere at spatial infinity is mapped to one point 
in field space, one can compactify the space to and 
consider mappings 5*^ — > SU{2), corresponding to pure 
gauge field configurations. 

One can also make other choice of gauge transforma- 
tion function 0(r) at spatial infinity (as long as the fields 
are pure gauge and constant in time there). In our case 
it is convenient to set {2n — 1)tt at r ^ oo. This 

choice, called "periodic instanton gauge" in this paper, 
in the original 4-dimensional theory corresponds to map- 
ping of the sphere at spatial infinity onto the equato- 
rial sphere of the SU (2) gauge group, parameterizing 
the pure gauge field configuration. This behavior of is 
equivalent to the requirement that the fields satisfy the 
following boundary conditions at ?■ = and r = oo, 



X|r->0 



(30) 



The conditions for the field at r ^ make the original 
field <i> regular at the origin. 

In this gauge no r-independent vacuum exists, but 
transition between vacua with n = and 7i = 1 is de- 
scribed in a very symmetric way. The behavior of the 
fields X and (j) for such transition is shown in Fig. |S1 In 
the original 4-dimcnsional model this topology changing 




r = 



FIG. 8: Topological transition in the SU{2) Higgs model: 
behavior of the fields (j) and X- Bold arrows show the change 
of the field as the radial coordinate increases from r = to 
r = oo. The configurations are shown: (a) at initial time, (b) 
in the middle of the process and (c) at final time. 



process corresponds to a transition where the fields wind 
over the lower hemisphere of SU{2) before the transition 
and over the upper hemisphere after the transition. 

The initial 0-boundary conditions in gauge theory are 
quite complicated. The basic reason is that there is a 
redundant field among the five fields oi, cj), (f>, X: Xj 
while the 6'-boundary conditions Ij9b|l are to be imposed 
on physical fields only. The analytic expressions for the 
modes /fc, gk in terms of the fields x, </>, ai are cum- 
bersome (see Refs. 0, [S^) and will not be presented 
here. It is simpler, and more precise in the lattice case, 
to perform this expansion numerically in the discretized 
version of the model. This expansion will be described 
in the following section. 

To complete the boundary value problem, one has to 
impose Gauss' law and the equation fixing the time in- 
dependent gauge invariance. Note that both of these 
equations are not full complex valued equations (unlike 
the 0-boundary conditions), otherwise the system would 
have been overdetermined. The point is that, the real- 
ity conditions at final time H9c|l forbid gauge transforma- 
tions with imaginary gauge functions and also guaran- 
tee that Gauss' law does not have imaginary part. So, 
only the real part of Gauss' law (|26a|) and equation fix- 
ing only real- valued gauge transformations must be used. 
Together with four 9 boundary conditions this gives the 
right number of boundary conditions for the system with 
five complex valued fields ai, a, f3, /x, v. The exact form 
of the gauge fixing condition will be given in section Hvl 
because it is again most conveniently expressed in lattice 
terms. 

One more complication of the problem is the invariance 
of the equations under translations along the real time. 
To solve the equations numerically this should be fixed in 
a controlled way, to make sure the contour winds around 
the branching points of the solution, and does not get 
too close to them. A method of removing this invariance 
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will also be described in Section ITVI 



IV. COMPUTATIONAL CHALLENGES 
A. Discretized action 

To obtain a self-consistent system of equations, the 
discretization of the equations (|26bP " (|26t|l should be done 
in a gauge invariant way. 

First, let us consider the discretized version of the ac- 
tion (|24|l . The spatial axis is discretized by introducing 
sites Ti, i = 0, . . . , TV, where tq = 0, r^v = L. The time 
grid consists of sites tj, j = —1, . . . , Nt + l. We arc work- 
ing in the ao = gauge, and omit the subscript in the 
spatial component of the gauge field, ai(r, i) = a{r,t). 



the Ari factors in the definition of a^-. The boundary 
conditions in the periodic instanton gauge, Eq. I|3U|) . are 



Xoj 

(f>Oj 



-MOj/2 



XNj 
(f>Nj 



i (31a) 

-i (31b) 

i (31c) 

-i (31d) 



for all j. In the boundary condition for (f>oj, the spa- 
tial derivative in Eq. H30|l was changed into a covariant 
one to preserve exact lattice gauge invariance. Thus the 
complex lattice field variables left are 



Field variables 



Xij: Xij 



and 



correspond to field 



Xij I Xij J Vij I Vij 

a. 



,^-1 .? = 



-1, ■ 
-1,. 



values on the space-time lattice sites, while de- 
fined on spatial links and temporal sites. We also absorb 



J 



S = {Sff + St + Si- + Sint) 

Nt N-1 

Sff = ^ vi.jWi^i{l - cos(a,j+i - g^j)) 

Nt N-1 



The discretized action reads 



(32) 



^ '"hAw3,^{X^,J + l - XiJ){X^,j + l - 

j=-l 1=1 ^ 

Nt N-1 , 

S'r = - XI ^%A^^A^'"'''X^+1,J - X^J){e~'"'''X^+l,J " X^J ) 

i=-l i=0 ^ 



Xij) + W4:^i{4>ij + i - 4>ij){4>i^j + i - 



+ «;i,,(e'''-/2^,+i,, - 0,,)(e-''^-/2^,+i,, - 0,,) 



Nt N-l , ^ 

j=-l 1=1 



T^iXtjXtj + ^)4>i]4>i] + 7^Xi]<t>'ij - 7;Xtj(l>l 



r 



where the weights are 



V2,j = hjAt^ 



Wl,i ''f+l/2/^''»+l/2 ^2.1 = l/Ari+i/2 W3.i 



rfAn 



with Atj^i/2 = tj+i — tj, Atj = (Atj+i + Atj)/2, and 
analogous expressions for Ar^; hj = 1 ior j = 0, Nt and 
1/2 for j = — 1 , A^i + 1 . Lattice field equations are derived 
from H32|l by variation over the lattice fields, after the 
boundary conditions Ij31|l have been used to exclude xoj j 
Xoj, <Po], 4>0] and XNj, Xnj, <l>Nj, (f>Nj from the lattice 
action. 



The action (|32|l is exactly invariant under time- 
independent lattice gauge transformations of the form 



a 

Xtj 
(t>ij 



^i+l — ^i J 



^i^i 



(33) 



This gauge freedom has to be fixed by boundary condi- 
tions. 
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B. Boundary term: normal modes 

To obtain lattice version of l)9b|l one notes that plane 
waves are no longer eigenfunctions of the Hamiltonian 
on the lattice. To find their analogue one brings the 
quadratic part of the action H32() . taken in the limit of 
continuous time, to the canonical form. We expand it 
near the space-independent vacuum H28(l 

X = -*-X X = i-^ (34) 

(j) = —1 + i<f> (j) ~ —1 — icj) 



(performing in the end a gauge transformation to the vac- 
uum H3U|) is straightforward). It is also useful to change 
to the notations 1)25(1 . 

X = i5 + i/3 x = a — i/3 

(j) = jl + iv (f) = fl — ii) 

In these terms, the quadratic part of the action ((32|l is 



5(2) 



' AT-l 



i=i 

2 / ~ ~ ^ ^ 

(a, + ai+i - ai) + - /?, 



I j=o 

{N-l 
2w2s 
i=0 

-I- ^ 2wi^i {ai/2 H- jli+i - jl{f H- {v^+i - 

N-1 N-1 

+ {Aw5,^ + W3,,)Pf + W3,» " 4ai/i, + 4/^2] -f 8A ^ 



N-1 



r 



(35) 



As seen from Eq. (|35|) , the variables ao , Po, P-o and i^o do where 
not have kinetic terms. Three of them are fixed by the 
boundary conditions at r = 0, 



ao = /3o = Mo = 0. 

The fourth one, vq, is determined from the field equation, 
which for this variable reads 

After the variables with 2 = have been excluded in the 
manner described above, the quadratic action takes the 
form 



Sij = -rSij — 
a J a J 



The symmetric matrix Sij is then diagonalized 

where Ok.t is an orthogonal matrix. Introducing yet an- 
other set of variables zj by the relations 

OkjUj = zk; yj ^ O'^j^zk , 
we finally bring the action to the diagonal canonical form. 



S 



where the real valued coefficients di and Sij are to be 

read oS^ from Eq. indices /, J label fields and space Therefore, vectors 
points, and (pi stands for the fields {cti, (3i, fliji^i^ai}. The 
change of variables 



dt 



1-2 1 2 2 
— Zr COtZt 

2 ^ 2 ^ ^ 



?(A')/ — OkI 



yi = dupi 

brings the kinetic term to the canonical form. 



S =^ I dt[^fj - ^yiSijyj ] , 



are normal modes in the lattice formulation of the theory, 
and should be used instead of the usual spherical waves. 
The corresponding frequencies are ujj^. 

The matrix Oxj and frequencies loi^ are found nu- 
merically. Since they depend only on the spatial lattice 
parameters (size and spacing) and coupling constant A, 
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and do not depend on the background vacuum field con- 
figuration, it is sufficient to perform this diagonahzation 
once for a given lattice. The first 4N — 3 eigenvectors 
^{K)i ^-nd eigenvalues luk correspond to physical modes, 
and the rest iV — 1 of them have w = and thus corre- 
spond to the gauge (unphysical) degrees of freedom. 



b. Zero modes part. The 0-boundary conditions H37() 
give only AN — 3 (complex) equations, while 5N — 4 
boundary conditions are required at the initial time. The 
"left over" N — 1 conditions correspond to TV — 1 gauge 
degrees of freedom in the model. As described in sec- 
tion for these iV — 1 equations one has to use the 
real part of Gauss' law 



C. Boundary conditions. 

a. 9 boundary conditions. To derive the lattice ver- 
sion of the boundary conditions, one takes the variation 
of the exponent for the total probability H15|l . which can 
be written in the following form 



iS{z) - iS{z') 



(36) 



where dots denote terms irrelevant in the current context 
and 7 = e~^. One has to vary the lattice version of l|3t)|) 
with respect to z/^-i (values of z at the first time slice) 
and set z' ~ z*. The variational equation reads 



SS 



5zi^-i 1 — 7^ 
which leads to 



(1 + 7^)^;/,-i + 



2t^/7 

1-72 



— — LOi tiezi^-i — - tj/lniz/._i = 0. (37) 



5zi^-i 1 + 7 



1-7 



where the derivatives of the action arc equal to the clas- 
sical momenta of the modes 



_5S_ 

Szi-i 



-ui._i(z7^o - -1) 



Here the index 7 = 1,..., AN — 3 labels physical degrees 
of freedom. One can go back to the original notations by 
means of the relation 

zi = ^.ijdjipj, 

where / = 1, . . . , 4iV - 3. 

Finally, we use the gauge transformation 



rii = TTCXp 



c(L-n) 



with c = 0.5 to transform the fields from the periodic 
instanton gauge l|Hn|l to the form (pj (of course, other 
gauge choices with il{0) = tt, ^l{L) = are possible). 
This transformation has the form 



(38) 



where ipj stands for the fields aj , f3j , fj,j , vj , Oj in periodic 
instanton gauge. The matrix gxj and vector (pY'^ can be 
easily read off from the expression for the lattice gauge 
transformation (|33|1 and definition (|34|l . 



Re sin(ai_i,_i - ai_i,o) - sin(aj__i - ai,o) 

+ iwsAxi-iXi.o - Xi-iXuo) 

1 — — 1 

+ 2*^4, -101,0 - 0i,-i0i,o) =0, (39) 

where i = l,...,A^— 1. One also makes use of equations 
that fix the remaining real gauge freedom. The latter 
equations arc 



Re = 



(40) 



for all L = AN — 2, ... , bN — A. These modes have zero 
frequency lol ~ Q and correspond to the unphysical de- 
grees of freedom which change under gauge transforma- 
tion, so Eq. H40|l fixes the residual gauge invariance with 
real gauge functions completely. Gauge transformations 
with imaginary gauge functions are forbidden by the re- 
ality conditions at final time. 

c. Final boundary conditions. It is straightforward 
to implement the reality conditions at final time (|9c|l . 
Supposing that the last two time grid points, Nt and 
A^t -f- 1 are on the real time axis, they are 



Ima^^TVt+i = Im/Jj^ATj+i = 



(41) 



For energies below the bifurcation line Ei{N) the time 
Nt can be chosen to coincide with the point C of the time 
contour (so there are only two lattice points Nt and Nt + 1 
on the whole CD part). For higher energies, though, the 
fields are not real along the most part of the real axis, 
so the part CD of the time contour has to be as long as 
possible (see Section IvTl. 

d. Fixing time translational invariance. One more 
complication is that, in the continuous formulation, the 
boundary value problem Q has an invariance under 
translations along real time (both field equations and 
boundary conditions are invariant under such a trans- 
lation). To define properly the boundary value problem, 
one has to fix the position of the solution in time. In 
the lattice version this invariance is violated by the dis- 
cretization and finite volume effects, but this violation 
does not enable one to control the position of the time 
contour relative to the branching points of the solution. 

The existence of this invariance means that one of the 
equations is redundant (if discretization and finite vol- 
ume effects arc discarded). Somewhat arbitrarily, we 
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take as redundant one of the real equations entering the 
0-boundary conditions (|^ 



arg/fc = arggfe , 



(42) 



for a specific mode. Provided the system hnearizes at 
initial time, this equation is indeed a consequence of the 
others. The reason is that reality conditions at final time 
imply that the (conserved) energy is real. Hence the 



J 



(linearized) energy H10() is real at initial time. Then one 
of the modes automatically obeys Eq. (|42|l provided all 
other modes obey the ^-boundary condition (|9b() . 

This suggests the following modification of the equa- 
tions. One of the equations (|9b|) is changed to 



I4l=e-''|5fe|, 
whose lattice version is (cf. Eq. ((23) 



(1-7^) 



SS 



^K-l 



2u) 



K 



Rc 



5S 



■ Im zx - 



Im ■ 



SS 



5zK, 



■ Rc _ 1 



(43) 



r 



Thus, instead of Eq. \A2\ one imposes another boundary 
condition, which is not invariant under time translations. 
The choice of the latter is a matter of convenience. We 
control the position of the solution in time by imposing 
the boundary condition that fixes the "center-of-mass" 
of the field x at the initial time to be equal to a given R: 



Re 



N-l 

E 



Wi.i{r.i - R){xi-iXi-i - 1)^ 



(44) 



This prescription works if the mode zk in Eq. (|43|l is rea- 
sonably occupied at the initial time, otherwise the equa- 
tion (|42|l . which is "thrown away", is nearly degenerate. 
Aside from this, the results of the calculations do not 
depend significantly on the mode chosen. 

The relative phase between fk and gk can be used to 
check the validity of the calculations. In the linear regime 
it should be equal to zero, so the actual value of this phase 
indicates how close the system is to the linear regime at 
the initial time. 

To summarize, the lattice boundary value problem con- 
sists of the field equations, obtained from action H32() for 
all inner lattice points {i = 1, . . . , N — 1, j = 0, . . . , Nt, a 
total of (A^— l)(A^t-|-l) equations), the final reality bound- 
ary conditions H41|) {N ~ 1 equations), the 6 boundary 
conditions (|37|l for all modes except one mode zk {N — 2 
equations), and a pair of real equations H43|) . (|44|) (one 
complex valued equation). This makes {N — l){Nt + 3) 
complex equations for the same number of variables. 



D. Search for solutions 

The equations to be solved make a set of discretized 
partial differential equations which change their signa- 
ture from hyperbolic on the Minkowskian parts of the 
time contour to elliptic on the Euclidean part. The prob- 
lem at hand is a boundary value problem which cannot 
be transformed into an initial value one. This means 



that the equations can be solved only globally, as a set 
of nonlinear equations at all r, t grid coordinates. 

To deal with the non-linear system of equations we em- 
ploy a multidimensional analog of the Newton-Raphson 
method which approaches the desired solution iteratively. 
At each iteration, the linearized equations in the back- 
ground of the current approximation are solved. The 
next approximation is obtained by adding the solution 
to the background; and the procedure is repeated. The 
advantage of the algorithm is that it does not require 
positive-definiteness of the matrix of second derivatives. 
It is, however, sensitive to zero modes. In the absence of 
zero modes, the algorithm converges quadratically; the 
accuracy of 10~^ is typically reached in 3-5 iterations. 
The convergence slows down in the presence of very soft 
modes, as typically happens near bifurcation points. 



E. Elimination algorithm 

The discrete version of the equations derived from H32|) 
dS 



IS 



d(pji 







(here ipjj ~ {a, a, /3, ^, i^}(<j, r^) and / runs from to 
5A^ — 4, = ... At). The Newton-Raphson iteration is 



d(pjidipj-i,K 



dipjidipjK 



+ Uj+i,K + TT— = (45) 



dtpjidtpj+i^K 



(all other second derivatives are zero) and before the next 
Newton-Raphson step the fields are changed according to 

(ri+l) (ji) 

Equations H45|) can be rewritten in matrix form 

D^^ ■ uj-i + Dj ■ uj + d'^^ ■ Uj+i +b^^0 (46) 
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where Uj and bj = are (57V— 4)-dimensional vectors, 

are (5-/V — 4) x (5A^— 4) matrices. I3y multiplying Eq. 
by we get 

u, = i?]-^ • u,-i + D^^ ■ uj+i + b, (47) 

with 1)]=^^ = -i)^! . bj = -DJ^ ■ bj. This system 

of linear equations was solved by the following version of 
"divide-and-conquer" elimination algorithm. Excluding 
Uj for some j gives 

+ {D^+\-b,+b,^,)_ 

+ {D\l\-b,+b,+,)^ 

Since the elimination of an equation changes only adja- 
cent equations, it is possible to eliminate all equations 
with odd J in parallel, and arrive to a system of the 
type (|47|l again, but with two times less variables and 
equations. This is the second level of elimination. After 
a series of eliminations we arrive at a system of only two 
equations for j = and j = Nt- 

uo = D^o~^ • it_i + D^i"^^ -UNt+bo (48a) 
"JVt = O'k^^ ■ uq + l)^^ • UNt+i + bMt (48b) 

where Z?'^-* and b have the values resulting from the elim- 
ination of all intermediate equations. Solving them to- 
gether with the boundary conditions* (also linearized), 
which involve u^i and ito for initial boundary condition 
and ujVt and u^t+i for the final one, we determine the 
corrections w-i, wq, ujVt and ujVt+i- Then it is straight- 
forward to reconstruct u at all intermediate points, using 
the equations H47|) for each elimination level. 



Unlike the field equations at the intermediate points 1451 , which 
are analytic, the boundary conditions involve complex conju- 
gation. So, equations I48i and boundary conditions are to be 
viewed as eight real matrix equations. All the elimination cal- 
culations (and reconstruction of field values afterwards) can be 
done, however, with complex algebra, which is two times more 
efficient. 



F. Solutions below the sphaleron energy 



The Newton-Raphson method requires a good initial 
approximation for the solution. This favors the follow- 
ing general strategy. We first find the periodic instanton 
solution, which corresponds to 9 = Q |l2l | and can be ob- 
tained via a minimization procedure. After the periodic 
instanton is found, we change the parameters T and 9 in 
small steps, using the solution from the previous step as 
a starting configuration. At each step we then calculate 
the energy E, number of particles N and the suppression 
exponent F{E, N) for the solution obtained. 

This procedure is illustrated in Fig. |51 where each dot 
represents one solution of the boundary value problem. 
Initial periodic instanton configurations correspond to 
the points on the upper left line in the figure. Start- 
ing from these points, the value of 9 was increased, and 
lines with constant values of T were obtained until the 
bifurcation line was met. Data obtained in this way make 
almost straight lines in the left part of Fig. |51 

The boundary value problem Q-lQlJ does not explic- 
itly refer to the topological properties. Hence, it is not 
guaranteed that its every solution describes a transition 
between topologically distinct vacua. This is not a prob- 
lem at = 0, because of the proper topological structure 
of the periodic instanton solutions. But at non-zero 9 one 
should check that the solution indeed has correct topol- 
ogy. 

The topological properties of a given solution are as- 
sociated with the behavior of the phases of the fields, 
see Fig. |S1 A very useful tool to control the properties 
of the solution is visualization of the field behavior. The 
visualization of a representative field configuration is pre- 
sented in Fig.^l It describes the field ^(''7^)1 with the 
phase of the field encoded in color. The Euclidean part of 
the time contour is inclined to make it distinct from the 
Minkowskian parts. In the initial state (left part of the 
surface) the field is close to its vacuum value, with exci- 
tation in the form of the incoming spherical wave moving 
towards to r = 0. The final state (right part of the sur- 
face) contains the outgoing wave. The phase of the field 
clearly behaves differently in the initial and final states. 
This confirms that the topological transition indeed has 
occurred (compare to upper three images in Fig.|HJ). Sev- 
eral other important properties of the solution may also 
be seen immediately. These are: the moment when the 
field goes through x = in the middle of the Euclidean 
evolution, which of course should happen with the field 
evolving between neighboring vacua; the wide outgoing 
wave, suggesting that a large number of low energy par- 
ticles is created after the transition; the small and rela- 
tively sharp incoming wave, meaning that higher energy 
modes are occupied and the number of particle in the 
incoming state is smaller. 
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FIG. 9; Search for solutions. Each point corresponds to one solution of the boundary value problem. The color of the points 
tracks the suppression exponent F{E, N). The almost vertical line is the line of bifurcations Ei{N), (cf. Fig.j^J. 



V. GOING OVER THE SPHALERON ENERGY 

The procedure described above works as it is for rel- 
atively low energies E < i?sph only. With growing en- 
ergy, the solutions on the CD part of the contour tend 
to stay for a long time close to the sphaleron. As the 
energy approaches some A^-dependent value Ei{N) this 
time tends to infinity, and if one continues to search for 
solutions to the boundary value problem (jSJ with reality 
condition imposed at finite positive time, the solutions 
above this energy have wrong topological properties, i.e. 
they end up in the same topological vacuum as the ini- 
tial one (see Fig. Ill|l . This situation is not specific to 
the SU{2) gauge model studied here, but appears quite 
generally in quantum mechanical tunneling with multiple 
degrees of freedom. It was observed also in the study of 
the false vacuum decay in scalar field theory Is^l and in 
quantum mechanics with two degrees of freedom 36, 33 ■ 
The phenomenon was studied in detail in Ref. j38| in the 
case of quantum mechanics of two degree of freedom, and 



a general method of dealing with this difficulty was pro- 
posed there and checked against the exact solution of the 
Schrodinger equation. We describe here its gauge field 
version. 

As suggested in Ref. [s^ , the line Ei (N) is the bifurca- 
tion line at which two types of solutions to the boundary 
value problem Q meet. These are (i) solutions which 
end up close to the same vacuum as the initial one and 
(ii) solutions that arrive at the sphaleron with excited 
positive modes (in the case of field theory these excita- 
tions fly away quickly in the form of spherical waves in 
the sphaleron background) . The former solutions are un- 
physical, while solutions of the latter form determine the 
tunneling exponent. For the interesting solutions of type 
(ii), the condition is satisfied only asymptotically, 
so it is very hard to find them numerically. A way out 
is to introduce a small regularization parameter into the 
equations of motion, which would not allow a solution to 
stay close to the sphaleron for infinite time. The final 
result is then obtained in the limit of zero regularization 
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FIG. 10: Visualization of the field x for a solution with = 1 
and E = 3.35. The color tracks the phase of the field. The 
part corresponding to the Euclidean evolution is inclined for 
visualization purposes. 




FIG. 11: Solution for T/2 = 2 and 61 = 3.35, without regu- 
larization. For this solution E/Esph = 1.04, Naw = 0.94, so 
that E > Ei{N). One observes that the topological proper- 
ties of the solution are wrong: it begins and ends in the same 
vacuum. 



parameter. 

To implement these ideas we start with the regularized 
expression for the cross section 

a,{E,N) = ^ |(/|e-2^^-5P£;PA,K>p , (49) 

where e is a small parameter and Tint is a functional 
proportional to the time the system spends in the inter- 
action region. In case of gauge-Higgs theory we use the 



functional 

Tint = Jdtjdr {4>{r)<P{r) - if . (50) 

The path integral for H49|) is no longer saturated by classi- 
cal solutions spending infinite time close to the sphaleron, 
while the original cross section <t{E, N) is obtained in the 
limit e ^ 0. The boundary value problem for (|49|l co- 
incides with the unregularized one, but with the action 
modified by adding an imaginary term of the form 

6S = ie jdt j dr {4>{r)(t){r) ~ i f . (51) 

The equations of motion (|26|l are modified accordingly. 

The functional (|50|l has several important features: i) 
it is gauge invariant; ii) it is large and positive on config- 
urations close to the sphaleron (where the fields stay for 
a long time away from their vacuum values); iii) it does 
not change the free dynamics in the linear region, since it 
docs not produce quadratic terms in the expansion of the 
action about vacuum (this is important for the boundary 
conditions to be unaffected). 

With this regularization one obtains results for all ener- 
gies E. The procedure is as follows. One introduces small 
but non-zero e at energies below Ei{N), then obtains so- 
lutions with proper topology for any energy Then one 
takes the limit e 0. Upon taking this limit, the con- 
figurations with E > El (N) stay for longer time close to 
the sphaleron, which means that in the limit of e ^ 
the solution tunnels "onto" the sphaleron. 

Moreover, at the boundary of the classically allowed 
region, the solutions to the regularized problem merge 
smoothly with the classical over-barrier topology chang- 
ing solutions, because the bifurcation on the boundary 
of the classically allowed domain is regularized exactly in 
the same way as the bifurcation at E = Ei{N). At the 
boundary of the classically allowed region F = by def- 
inition, so the regularized version of this functional i^e is 
proportional to e. This means that T and 9 are also pro- 
portional to e there, and as the regularization is turned 
off, e ^ 0, both T and 9 disappear, leading to purely real 
classical boundary value problem in real time. 

There is one more complication in the SU (2) field the- 
ory, which is relevant to this procedure. For energies 
E < El (N) the amount of time T spent on the Euclidean 
part of the contour is a growing function of energy (as 
opposed to the situation in two dimensional quantum 
mechanics H^), while it is zero at the boundary of the 
classically allowed region Eo{N). This means that T as 
function of energy at fixed 9 has a maximum somewhere 
in between Ei{N) and Eq{N) (the calculations show that 
the maximum is actually at Ei{N), see Fig. 113(1 . For 
T close to this maximum value, the Newton-Raphson 
method fails because of the presence of two nearby so- 
lutions with equal values of T. This new bifurcation 
is absent, if one searches for solutions with fixed energy 
E, instead of fixed T. To formulate the boundary value 
problem with fixed E instead of fixed T, one simply sets 
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the contour height T' to some conveniently chosen value, 
and leaves T as a free variable. This leads to a trivial 
modification of the initial boundary conditions, 

,4 = c-«--^(^-^')g, . (52) 

An additional equation is then required to find one extra 
undetermined variable T . This is the equation involving 
energy of the solution, Eq. Hl()|l . With this modifica- 
tion of the procedure, the "bifurcation" corresponding to 
the maximum of T disappears. Note however, that this 
method had to be applied with great care — on realistic 
grids it is hard to achieve fully linear regime in the initial 
state, therefore the difference between T and T' must be 
small, as the dependence of the fields on imaginary time 
is exponential. In our calculations, the version of the 
boundary value problem with fixed E (T' different from 
T) was used only to cross the maximum of T . 

In our numerical calculations we introduced small non- 
zero e as the energy of the solutions approached Ei (N) . 
Simultaneously, the modification H52|) was used to get 
past the maximum of T at Ei{N). In Fig.|51the solutions 
obtained with this modification are represented by points 
on the line that crosses the bifurcation line Ei{N). At 
higher energies the modification (|52|l is no longer needed, 
and only regularization with non-zero e was used. In 
Fig. |51 solutions to the regularized problem in region A. II 
correspond to points on the curved lines in the right part 
of the plot (lines of constant T). The line with the highest 
energy has zero suppression exponent and corresponds 
to the boundary of the classically allowed region. To 
connect lines of constant T, we obtained a set of solutions 
represented by the irregular line in right part of Fig. |3| 

VI. NUMERICAL RESULTS 

There are several factors affecting the choice of the lat- 
tice size and shape. The physical spatial size of the lattice 
L is chosen large enough to make comfortable room for 
the sphaleron. More importantly, L determines how close 
to the linear regime the system is in the initial state: the 
farther away from the origin, the smaller becomes the 
amplitude of incoming spherical waves. After L is cho- 
sen, the lengths of the parts AB and CD of the time 
contour are determined completely; the length of the AB 
part Tab is slightly smaller than L, so that the incoming 
wave does not reach the spatial boundary r = L. The 
length of the CD part Tcd is zero for energies below the 
bifurcation energy E < Ei{N). For higher energies, Tcd 
is adjusted to be long enough, so that the solution gets 
close to the vacuum configuration, and the regulariza- 
tion (|51|l docs not contribute significantly to the equa- 
tions of motion at the final moment of time when the 
reality boundary conditions are imposed. 

The lattice spacing Ar constrains the precision of the 
discretization in two different ways. First, it is chosen to 
be substantially smaller than the size of the instanton- 
like part of the configuration, i.e. the characteristic scale 



of nonlinear dynamics that occurs near r = during 
the topological transition itself. Second, Ar controls the 
energy of the hardest mode in the initial state, thus lim- 
iting the lowest particle number N that can be reached 
for given energy E. The time spacing At is chosen to be 
smaller than Ar to guarantee stability of the numerical 
procedure. 

The amount of computer memory required for a lattice 
of spatial size Nr and time length Nt is approximately 
2 X Nt {5Nr)^ X 16 bytes (see Section llVEj) . while the CPU 
time of one Newton-Raphson iteration scales roughly as 
Nt{5Nr)^. It was noted in Section ITV El that the algo- 
rithm is suitable for parallel execution, so one divides 
this time by the number of processors available for the 
calculations^. Overall, the most strongly constrained is 
the spatial size A^^: a two times larger spatial grid means 
eight times longer processor time. 

The main results we present in this paper were ob- 
tained on a grid with spatial size L = 8 (i.e. L = 
8/y/2Mw) and number of spatial grid points Nr = 90. 
The length of the initial Minkowskian part of the con- 
tour Tab was equal to 6. The number of time grid 
points Nt on the part AB of the contour was equal to 
200, while on the Euclidean part BC it was equal to 150. 
The number of points on CD part varied from 2 for ener- 
gies E < Ei{N) to about 400 for higher energies (when 
the e-rcgularization was used). On the largest grids the 
amount of memory used was 4Gb, and it took 3 min- 
utes for one Newton-Raphson iteration on a 16 proces- 
sor IBM-RS/6000 supercomputer, or about 15 minutes 
for one full solution. 

We obtained the results for the suppression factor in 
the region of E and TV shown in Fig. For the lattice 
parameters we used, this region is limited mainly by the 
effects of nonlinearity at the initial time, preventing us 
from reaching smaller particle numbers. When energy 
and particle number are small simultaneously (bottom- 
left part of the plot), effects of the spatial discretization 
(finite Ar) are also important. 

To check the discretization effects, a limited set of 
calculations was performed on smaller grids. The re- 
sults presented here coincide with results obtained with 
Nr = 64 with precision better than 1% (except for very 
small energies). With Nr = 45, on the other hand, the 
results coincide only for sufficiently large initial particle 
numbers, exactly as one would expect. 

The linearization of the system in the initial state can 
be checked by evaluating the time dependence of the lin- 
ear energy (|10|l and particle number (|ll|l on the part AB 
of the contour. For linearized system, these should be 
independent of time. For a typical configuration this test 
is shown in Fig. 1121 The linear energy coincides with 



The parallelization algorithm is effective only if A''procossors < 
V Nt , so the shortest possible wall clock time in an ideal situation 
is proportional to y/ Nt{5Nr)'^ 
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FIG. 12: Particle number evolution (upper plot) and linear 
energy evolution (lower plot) for the configuration with N = 1 
and E = 3.35. The exact (full nonlinear) energy is also plotted 
for reference (straight dotted line). 



the exact one in the initial state with precision of or- 
der of 1% or better, which confirms that the solution is 
quite close to the linear regime. Another test of linear- 
ity is the amount of the violation of the initial boundary 
condition H42|) . which is discarded to impose the time 
translation invariance fixing relation (Section lIVp . This 
amount grows towards smaller TV, and apparently this 
is one of the effects preventing us from going to lower N 
with the current spatial lattice size L = 8. Larger lattices 
are needed to achieve better linearization on the initial 
part of the time contour and thus reach smaller particle 
numbers. 

We made additional checks of the precision of the 
numerical calculations, including conservation of energy 
and the inverse Legendre transform 



dF 
^ ' dN 
dF 

^ - -^'-^ BE 



(53) 
(54) 



N 



These checks are satisfied with precision better than 
10"'^. This means that the precision of the final results 
is determined mostly by the quality of the linearization 
in the initial state (of the order of 1%). 

Lines of constant T and constant 9 are shown in 
Figs. 1131 and 1141 One observes from Fig. El that 9 grows 
as N decreases, as expected, and 9 is equal to zero on the 
periodic instanton line and on the boundary of the classi- 
cally allowed region Eq [N) . The lines of constant T show 
that T also equals to zero at the boundary Eq{N), and 
reaches a maximum for given TV (and for given 9 also) at 
the bifurcation line Ei{N). Close to this line we made 
use of the modification of the boundary value problem 
described at the end of Section El 

Representative solutions are shown in Fig. 1151 They 
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FIG. 13: Lines of constant T. 
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correspond to deep tunneling regime (E < Ei(N)), tun- 
neling onto the sphaleron {E > Ei{N)) and classical 
overbarrier transition at Eo{N), all for A'' = 1. One 
can see from the color patterns that the field indeed un- 
dergoes the topology changing transition of the form il- 
lustrated in Fig. |S1 The incoming wave is present in the 
left part of the pictures, becoming sharper and sharper 
for higher energy (the particle number is the same for 
all plots). In the first picture the topological transition 
is seen on the Euclidean part of the contour. In the 
second and third pictures, a sphaleron like configuration 
is visible on the right, with "extra" waves (excitations 
about the sphaleron) flying away, while the sphaleron it- 
self starts to decay quite close to the right end of the plot 
(with the regularization parameter e tending to zero, the 
moment of sphaleron decay moves towards larger times) . 
At large times, the wave reflected from the boundary 
r = L appears due to the Dirichlet boundary condi- 
tions (|31|) imposed at r = L. This wave does not alter 
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FIG. 15: Surfaces describing the x field for solutions with 
TV = 1 and -B = 3.35 (upper picture), E = 4.48 (middle), 
E = 5.22 (lower picture). The first surface corresponds to 
deep underbarrier tunneling, and the last one corresponds to 
nearly classical over-barrier transition. 



the results, as it occurs in the hnear regime^. 

The fact that for E > Ei(N) the solution after tun- 
neling has the form of the sphaleron plus spherical ex- 
citations in its background is illustrated by plotting the 
spatial energy density at different times after tunneling. 
In Fig. 1161 the energy density distribution is shown for 
the middle solution of Fig.^l As the time increases, the 
bump on the right (spherical wave) moves towards larger 
r, while the energy density profile approaches that of the 
sphaleron. 



To get rid of the reflected wave a much larger spatial grid would 
be needed. 



0.2 




r 



FIG. 16: Energy density for several values of Ket (indi- 
cated by numbers near graphs) for the middle configuration 
of Fig. iT^(iV = 1, E ^ 3.35). The energy density for the 
sphaleron solution is shown for comparison by solid line. 



It is also instructive to see that with the number of 
incoming particles decreasing, the occupied modes have 
higher frequencies. This is demonstrated in Fig. 1171 for 
energy slightly smaller than the sphaleron energy. 

Finally, let us discuss the extrapolation of the results to 
zero number of particles, which we performed to obtain 
predictions for the suppression exponent Fhg (E) of the 
two-particle cross section. 

Two ways of obtaining the lower bounds on Fhg{E) 
were explained in Section |n] One is to continue F{N) 
at each energy linearly in TV to = (this is justified 
by recalling that dF/dN oc —6 increases as A ^ 0), 
while the other is to continue lines of constant F (see 
Fig. 121) linearly to A = (this gives a lower bound since 
the lines of constant F have positive curvature). Both 
these extrapolations are straightforward to make, insofar 
as the required derivatives of F{E, N) are given for each 
configuration by the values of T and 9 through the rela- 
tions H53|l and (|54|) . In this way we obtained the bound 
shown in Fig. |S1 

We now elaborate on our estimate of the function 
Fhg(E) itself, Fig El Perturbative calculation at low en- 
ergies |43 shows that while the exponent F{E, N) has 
singular behavior of the form Nlog{N), the function 
T{E, N) is regular in A^ and close to a linear function. 
At large A^ the numerical data also demonstrate that the 
behavior of T is close to linear. This is shown in Fig. 1181 
For E > Esph numerical results suggest that T{N) is 
almost constant at small A^. These properties justify a 
linear extrapolation of T{N) to N = with energy kept 
constant. After obtaining T{E) at zero particle number, 
Fhg is readily found by integrating Eq. H54|) starting 
from the instanton value Fhg{E = 0) = 1. The result- 
ing estimate is presented in Fig. O solid line. 
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FIG. 17: Distribution of particle number rife at initial time 
over modes for Eaw /Mw ~ 3.54 and different A^. at are 
numbers of particles in each mode in units of 1/aw for four 
different types of modes (see Ref. IZH for definitions). Mode 
ai is the Higgs boson mode, while 02,3,4 are gauge boson 
modes (02,3 are transverse, 04 is radial). On the horizontal 
axis is the mode number for a lattice with spatial size r — 8. 



VII. CONCLUSIONS 

Our study shows that the semiclassical procedure cou- 
pled to suitable computational techniques is capable of 
producing quantitative results for semi-inclusive, weakly 
coupled non-perturbative processes, and, in particular, 
for particle collisions. 

In this paper wc applied this technique to study the 
suppression factor for topology changing transitions, and 
accompanying baryon and Icpton number violation, in 
the SU(2) sector of the electroweak theory up to ener- 
gies well above the sphaleron energy. We imposed spa- 
tial spherical symmetry, so our results are valid, strictly 
speaking, for s-wave scattering. 



Our results show that the known analytic expression 
for the suppression exponent, which contains three terms 
of low-energy expansion, works well up to the sphaleron 
energy, but underestimates the suppression at higher en- 
ergies. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 




FIG. 18; Left: T{N)/2 for different energies, labeled by values 
of Eaw /Mw , the points are data from numerical calculation 
and lines are extrapolations. Right: T{E)/2 extrapolated to 
zero particle number. 



By numerical analysis we have found that baryon and 
Icpton number violation, accompanying topology chang- 
ing s-wave particle collisions in the electroweak theory, 
remains highly suppressed up to energies of at least 
^ 250 TcV (and likely much higher). 
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